theta01 = -1.1e-16
theta11 = 0.88
theta21 = -0.05

uy = 3.4e5
sy = 1.23e5
ux1 = 2000
sx1 = 795
ux2 = 3.17
sx2 = 0.76

theta0 = (theta01 - theta11 * ux1 / sx1 - theta21 * ux2 / sx2) * sy + uy
print(theta0)
theta1 = sy * theta11 / sx1
print(theta1)
theta2 = sy * theta21 / sx2
print(theta2)

print('th1', theta11 * ux1 / sx1)
print('th2', theta21 * ux2 / sx2)
